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Two-level systems and mass deficit in quantum solids 
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We study a disordered quantum solid incorporating two-level systems in which a group of atoms (or a 
single atom) can experience coherent tunnelling between two different positions and demonstrate that an effec- 
tive mass deficit induced by the presence of such objects can manifest itself only at relatively high frequencies 
and should vanish in the low- frequency limit. The crossover to the regime which can be associated with the 
appearance of an effective mass deficit has been observed in recent torsional oscillator experiments. 
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1. Introduction. After a number of torsional oscil- 
lator experiments [1 - 3] demonstrated that at low tem- 
peratures solid 4 He behaves itself as if it contained a 
superfluid component which remains at rest when an 
external force is applied to the sample, an idea was put 
forward that these properties can be explained by the 
presence of two-level systems (TLSs) in which a group 
of atoms (or a single atom) experiences classical [4] or 
quantum ^ [5] tunnelling between two localized posi- 
tions shifted with respect to each other. To some ex- 
tent this conjecture is confirmed by the dependence of 
the experimentally observed mass deficit on the degree 
of disorder in the solid [2]. 

However, the analysis of Refs. [3 [6] has led to a 
rather paradoxical conclusion that the effective mass 
deficit induced by TLSs should remain non-vanishing 
even in the adiabatic limit including the really station- 
ary situation of the thermodynamic equilibrium, when 
experimentally the system behaves as if p s = [7| . The 
origin of this conclusion can be traced to the assump- 
tion that the bulk velocity of a solid is a classical variable 
commuting with the Hamiltonian of a TLS. 

In the present note we propose an alternative ap- 
proach which takes into account the operator nature of 
different variables in a more consistent way. This al- 
lows us to show that a quantum solid with incorporated 
TLSs can demonstrate the presence of a pronounced 
mass deficit only when an external force changes fast 
enough, but not in a stationary situation or at very low 
frequencies. For comparison, we also discuss the situa- 
tion when the tunnelling inside the TLSs is incoherent. 

2. A solid with a quantum two-level system. 
Consider a solid of total mass M incorporating a TLS 
in which a group of atoms (or a single atom) of mass 
m can tunnel between two localized positions shifted by 
vector a with respect to the rest of the solid. Below 



we use the term "solid frame" to designate the part of 
the solid which does not participate in the process of 
tunnelling and the term "system" to describe the whole 
system consisting of the solid frame and the TLS. 

In the absence of any interaction between the TLS 
and other internal degrees of freedom (for example, 
phonons) the quantum-mechanical Hamiltonian of such 
a system can be written as 



1 9 

H = — P 2 - ea 3 + Jen 
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where P = —ihd/dH is the operator of total momen- 
tum of the system (conjugate to the center of mass po- 
sition R), the second term describes the difference in 
energy (given by 2e) between the two localized states of 
the TLS and the third one the process of quantum tun- 
nelling (with amplitude J) between these two states, b\ 
and (73 being the Pauli matrices. 

Naturally, the velocity of the center of mass of the 
system, 

is determined just by its total momentum P and is in- 
sensitive to its internal life (a particular state of the 
TLS). In the presence of an external force F(t) applied 
to the system as a whole [which corresponds to replac- 
ing H by H = H - RF(t)], the time evolution of (V) 
is entirely determined by relation (d/dt)~V — F/M. 

However, in a number of experimental situations an 
external force (for example, of mechanical origin) is ap- 
plied not to the center of mass of the system, but to the 
solid frame, and one is interested in the relation between 
this force f(t) and the velocity of the solid frame 2 ', v(t), 
which is not obliged to have exactly the same form as 
the relation between V and F. The coordinate describ- 
ing the position of the solid frame, r, can be introduced 
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2 'Andreev calls the same variable "the solid bulk velocity" |51I5). 
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by rewriting the definition of the center of mass position 
R as 

MR = (M - m)r + tox = Mr + m(x - r) , (3) 

where x is the position of the center of mass of the atoms 
(atom) forming the TLS. After replacing u = x — r, the 
displacement of the TLS with respect to the solid frame, 
by — iao"3, Eq. can be rewritten as 



r = R 



2M 



acr 3 . 



(4) 



The operator of the solid frame velocity v = (d/dt)r is 
then given by 



H ,r 



M 



mj 



-acr 2 



(5) 



Naturally, this equation can be also rewritten as 

, , to Ja „ 

P = Mv —a 2 , 



(6) 



which corresponds to splitting the total momentum of 
the system P into the two terms related respectively 
with the solid frame and the TLS. However it is im- 
portant that in contrast to the center of mass velocity 
V = P/M, which in the absence of external force com- 
mutes with the Hamiltonian, the velocity of the solid 
frame v is an operator which does not commute with 
H . 

In the presence of an external force f(i) applied to 
the solid frame the Hamiltonian of the system acquires 
form 



H = H - rf(t) = -Rf (t) - h a (t)a D 



where 



h a {t) 



-J, 0, 



2M 



af(t) 
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plays the role of the effective magnetic field acting on 
the spin 1/2 which can be associated with the TLS and 
subscript a — 1,2,3 denotes the components of h. 
It follows from Eq. $5$ that 



dt 



f to Ja d 

M + Mlldi^ 2 ' 



(9) 



which suggests that the coefficient of proportionality be- 
tween (d/dt)(v) and f can be different from 1/M. To 
find a more explicit form of a TLS-induced correction 
to the solid frame's equation of motion one needs to 
express the time derivative of {a 2) in terms of f(t). 

3. Adiabatic regime. If external force i(t) does 
not depend on time or changes very slowly one needs to 
take into account the relaxation processes which force 



vector {a a ) to remain always parallel to h a , from where 
(er 2 ) = 0. Therefore, when f(t) changes sufficiently 
slowly the processes inside the TLS make no correc- 
tions to P = Mv. Thus, in the adiabatic regime the 
presence of the TLS cannot lead to the appearance of 
any difference between the real mass of the system and 
its effective mass observed in experiments involving the 
application of mechanical force. 

The finite expression for the mass deficit in the sta- 
tionary regime derived by Andreev [5] in the easiest way 
can be reproduced by calculating the average of the TLS 
contribution to the total momentum P, that is of the 
second term in Eq. ©, with the help of the Hamilto- 
nian, 



h M 2 toJ . 
H=-v -— (av)a 2 



£(73 + J&x , (10) 



which in Ref. [5] is obtained by applying to 
Htls — ~£&3 + J&i the Galilean transformation 
from the reference frame in which the solid frame is 
at rest to the reference frame moving with velocity v. 
Naturally, the application of this procedure implies 
that solid frame velocity v can be treated as a classical 
variable. 

One can easily check that Eq. (fT0|) differs from Eq. 
([1]) with P replaced by Mv — (raJa/S)^ only by a triv- 
ial constant term which additionally tends to zero in 
the thermodynamic limit, m/M — > 0. Thus the basic 
difference between the two approaches consists only in 
choosing whether P or v is a classical variable commut- 
ing with the Hamiltonian. The first option (adopted 
here) leads to (02) — and (P) = Mv, whereas the 
second one produces for the average of the TLS contri- 
bution to the total momentum the expression [5] 



TLS/ 



toJ 



tanh[E(v) /T] 

W) 



a(av) (11) 



with 



E(v) = \Je 2 + J 2 + [mJ/hf (av) 2 , (12) 



which does not respect the Galilean invariance being 
nonlinear in v. 

4. Finite-frequency response. When interaction 
of the TLS with other degrees of freedom (heat bath) is 
weak, the basic form of the linear response of the TLS 
to the application of the external force with a finite fre- 
quency, f(i) oc cos(wt), can be found by constructing 
the periodic solution of the equations describing the free 
evolution of the TLS El , 



dt 



[H,a a 



(13) 
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After diagonalizing Htls = —£03 + J&i and construct- 
ing the corresponding finite-temperature density matrix 
p = cxp ^— H^ls/T^J one obtains that in the absence of 
an external force 



do) _ 



-£,0, J)tanh| (14) 



with £ = \/e 2 + J 2 . 

Solution of the equations for (a a ) obtained by the 
linearization in the vicinity of (a a ) = (<7 Q )(°) gives 



dt 



af(<) 



Muj 2 -n 2 h 



>(0) 



(15) 



with fi = Substitution of Eq. (JUD into Eq. (15 

and then into Eq. ([9]) transforms the latter into 



d 

~dl { 



M 



aia,j 



CO' 



n 2 m 2 



m, (i6) 



where subscripts i and j denote the components of vec- 
tors in the real space and 



m2j \ ^ E 



(17) 



When one additionally takes into account the pro- 
cesses of transverse relaxation induced by the inter- 
action with the heat bath 3 ), the singularity at uj = 
57 is smeared out and the TLS-induced correction to 
(d/dt)(vi) acquires also a dissipative contribution, pro- 
portional not to cos(wi) but to sin(wi). If for simplicity 
one assumes that the relaxation can be characterized by 
the same relaxation time r (with Qt 3> 1) for both di- 
rections perpendicular to (a a )(°\ the form of the result 
corresponds to the replacement of uj 2 — fl 2 in the denom- 
inator in Eq. (fT6|) by [w + i/r) 2 — Q 2 . The dissipative 
contribution is the dominant one when t\u> — f2| -C 1. 

When the system contains many TLSs with ran- 
dom orientations, Eq. (j 16(1 can be replaced by 
(d/dt){v) = i(t)/M cS (u) with 
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M cS (uj) M 3M 2 ^ (uj + i/r n ) 2 - ft 



(18) 

where subscript n numbers different two-level systems 
and we have assumed that vectors a„ have an isotropic 
distribution. As usual with such a notation, the real and 
imaginary parts of 1/M e s(u)) describe the amplitudes of 
terms in (d/dt)(v) which are proportional respectively 
to cos(wi) and sin(u;i) when an external force applied 
to the solid frame is proportional to cos(wi). 

3 'The interaction with the heat bath can also lead to the renor- 
malization of J [S]. 



It follows from the structure of Eq. I|18p that when 
u> is much larger than all Q n the form of M e s(u>) indeed 
corresponds to the presence of frequency-independent 
mass deficit (as it was suggested by Andreev jS]). With 
the decrease in frequency the value of mass deficit de- 
creases and passes through zero at frequencies at which 
the dissipative contribution is most prominent. At 
uj — * both parts (dissipative and nondissipative) of 
the TLS-induced correction to l/M tend to zero. 

5. The case of incoherent tunnelling. Con- 
sider now the case when the tunnelling process inside 
a two- level system is incoherent. In such a situation 
the displacement of the TLS with respect to the solid 
frame u = x — r induced by a time-dependent external 
force f(t) cx cos(wt) coupled to the solid frame's position 
r = R (m/M)u acquires a very simple form [5], 



u{uj) = 



1 



a(af) 



-iruj + 1 M 4Tcosh 2 (e/T) 



(19) 



where relaxation time r is inversely proportional to the 
tunnelling rate, whereas e retains the same meaning as 
above. 

Substitution of Eq. (fH))) into the classical analog of 
Eq. (0, namely 



dt 



1 

M 



m d du 
Mdi~dt 



(20) 



then leads to the following expression for the frequency- 
dependent coefficient of proportionality in the relation 
(d/dt)v = f/Meff: 



1 



1 



£ 



i 



2 ? 



-ir n uj + lcosh\e n /T) 



M eS (uj) M Y1TM 2 

(21) 

The form of Eq. (f2Tj) demonstrates that in the case of 
incoherent tunnelling the TLSs-induced contribution to 
the effective mass can never be described in the form 
of frequency independent mass deficit and at high fre- 
quencies is of the dissipative nature. 

6. Conclusion. In the present note we have ana- 
lyzed how the presence of TLSs influences the dynamic 
properties of a solid and have demonstrated that in 
order to observe a reduction of the effective mass of 
the sample (analogous to that in superfluids) the solid 
should incorporate quantum TLSs and the frequency 
at which the external force is applied should be high 
enough in comparison with resonance times of the TLSs. 
In the case of incoherent tunnelling the regime in which 
the contribution from the TLSs can be described as 
frequency-independent effective mass deficit is absent. 
Since TLSs are local objects, analogous conclusions are 
applicable also to the effective moment of inertia 7 e ff(w), 
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which is a relevant quantity in torsional oscillator ex- 
periments. Experimentally, the crossover to the regime 
which can be associated with the appearance of an ef- 
fective mass deficit has been demonstrated in a number 
of works starting from those of Kim and Chan p.] . The 
detailed comparison of the temperature dependences of 
both components of the response can be found in Ref. 
[TO] , which also contains a comprehensive list of refer- 
ences to other experimental and theoretical works. 

An essential feature of our results is that in both 
regimes (of coherent and incoherent tunnelling) the 
TLS-induced contribution vanishes in the limit of 
lu —>■ Q. Another common feature of the two cases is 
that the presence of the TLSs makes an additive con- 
tribution to 1/A/ C ff or 1/ieff rather than to M e ff or I e ff. 
In that respect the situation is quite analogous to that 
in superconducting vortex glasses where TLSs make a 
positive contribution to the inverse superfluid density 
(in other terms, specific inductance) rather then a neg- 
ative contribution to the superfluid density itself [9] [TTJ . 

It seems worthwhile to mention that the phenomeno- 
logical approach [12] used in Ref. [10] for analyzing the 
experimental data is based on conjectures which are in 
contradiction with both these properties. Namely, the 
authors of Ref. [12] have assumed that the internal de- 
grees of freedom of solid 4 He make a contribution to its 
back action which in terms of the frequency-dependent 
effective moment of inertia of 4 He sample, I e ff(u>), can 
be written as an additive correction to the frequency- 
independent moment of inertia of the solid frame, 

9a 



/off M = i + 



(22) 



with /3 < 1. Moreover, at low frequencies this contribu- 
tion does not tend to zero but behaves itself like a neg- 
ative correction to the stiffness constant of the torsional 
oscillator. If it were really so, the solid by itself (outside 
of the torsional oscillator) would be unstable. A more 
logical assumption on the form of I c g(iu) in a glassy 
system with incoherent tunnelling inside the TLSs [con- 
sistent with the form of Eq. (j2"T]) ] would be 



9^ 



(1 



(23) 



iTU))P 

The main argument of Ref. 10] in favor of the su- 
perglass state consists in the impossibility to reconcile 
experimental data with dependence ([22]) - the observed 
frequency shift is too large in comparison with the maxi- 
mum in dissipation (characterized by the inverse quality 



factor), the same being true also for dependence (j2"3"j) . 
However, the dispersion of the parameters of the TLSs 
can be taken into account by the replacement of /3 = 1 
by (3 < 1 only for a particular form of the distribution 
and in other cases may lead to different dependences, 
especially if the appearance of the frequency shift is in- 
duced by the crossover between the regimes of incoher- 
ent and coherent tunnelling inside the TLSs. The tem- 
perature dependence of g, the amplitude of the TLS- 
induced term, also may change the relation between the 
frequency shift and the minimal quality factor. To clar- 
ify the situation, the experimental investigations of the 
temperature dependence of the solid 4 He dynamic re- 
sponse have to be complemented by more systematic 
studies of its frequency dependence. 
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